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ABSTRACT 


This thesis provides the student of combat modelling 
with a computer program for a relatively simple combat 
model that can be used in a classroom environment for study 
and analysis. The model is an aggregated, force-on-force 
ground-combat model that uses Lanchester's aimed-fire 
equations for casualty assessment. The original version 
of the model was developed in 1979 in a previous thesis, 
and shortfalls in the original model have been overcome, 
along with the addition of several enrichments to provide 
added user flexibility. A user's manual is provided to 
facilitate user TES to the model from a permanent disk 


шаі спе W. R. Church Computer Center. 
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л ЇЕ 2 


As the art of combat modelling becomes more advanced, 
combat modelers are continually building more and more 
complicated models. To the beginning modeler, the ability 
to understand how these models operate is difficult, if 
not impossible. As a student of combat modelling, 1 have 
sensed a need for a relatively simple model that could be 
easily studied in a classroom environment. It was with 
this in mind thst this project was initiated. 

At present, there seems to be no simple combat model 
available that demonstrates the basics of model building 
to the beginning student. In 1979, J. Smoler (a student 
at the Naval Postgraduate School) attempted to build such 
a model for his thesis research (Smoler,1979), His model 
was a deterministic, force-on-force computer model that 
used Lanchester's aimed-fire equations for casualty ass- 
essment. The general scheme of his model is shown in 
Ensure 1. 

Although Smoler's model was a bold attempt at a simp- 
listic combat model, it did have some problem А 
warranted investigation. Some of these problems were 
first discovered during a class project in a combat models 
class (OA 4655) at the Naval Postgraduate School. It is 
the purpose of this thesis to study Smoler's model and to 


undertake an enrichment program that will make the model 
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FIGURE 1. GENERAL SCHEME 





more realistic and flexible, while maintaining it's trans- 
parency and simplicity. 

Once enhanced, the model would then be made available 
for use аз an instructional tool for combat modelling 
classes. This will be accomplished through the develop- 
ment of a user's manual to facilitate the use of the model 
by students by placing the model on a permanent disk in 
the We Re Church Computer Center where it will be easily 
accessed by any desired usere 

The remainder of this paper will discuss, in detail, 
the problems found with Smoler's model and the methods 
used to solve these problems. It will also outline some 
new features that have been incorporated into the model 
to allow user flexibility. Finally, the User's Manual 
(Appendix A) will provide all the required information 


to enable even the novice user to utilize the model. 





ІІ. THE ORIGINAL MODEL 


A. GENERAL DESCRIPTION 

Smoler's original model of land combat is a deter- 
ministic model that plays combat between two homogeneous 
forces, a blue force and a red force. The blue force is 
comprised of three subunits in a static defense, with each 
subunit armed with three TOW antitank missile systems. 
The red force is composed of three subunits of three tanks 
each, attacking on pre-planned routes. The battle takes 
place on the 10 x 10 Km Fulda Box that has been developed 
and used in the STAR simulation model (Wallace,1978), Since 
eee Major thrust of this paper is to alter the original 
model to a form that will be easily used for classroom 
instruction, a brief discussion of the major components 
of the original model will be presented, including problem 


areas that have motivated model changes. 


В. ATTRITION PROCESS 

The attrition process in the original model utilizes 
Lanchester "aimed-fire" equations with variable attrition 
coefficients. The Lanchester equations used are fairly 
simple and will be discussed later. However, the calcu- 
lation of the attrition coefficients is of more immediate 
interest. Smoler used the Bonder-Farrell model to calcu- 


mabe the coefficients, A. 


ij the rate at which one firer 





of unit i kills unit j targets. These А + 5 are computed 


eceording to: 


A... = МЕ.) 


Hg 
where EU 4) is the expected time for one firer of unit i 
ЕЗ 111 опе target of unit j. The B(T, j) is computed 


using the Bonder-Farrell formula: 


ЕСТ. 4) = es + Ч h 


Ct t.)/P(h Im) x ((1-P(h|In))/P(KIE) + 


тэ. 1. t.)/PCKIE) + 


P(h|h) - P) 


where 
t = time to acquire a target 
Ч eco Lice first round following acquisition 
t ео tire following a hit 
$ = time to fire following a miss 
Us = time "ome ent of a round 


шэг 1202 Of a first round hit 


ИИО роо аа Лау ое а nit following a hit 
БЕС probability of a hit following a miss 
POIS probabi lites of a kill given a hit 


Шие Тогпціа holas гот the conditions that the hit 
probability of any round depends only on the result of the 
previous round and no accumulated damage is considered. 
smoler also assumed that for the TOU weapon system, 


P(KIH) = 1.0 and P(hlm) = P(h|h) = P, thus reducing the 


10 





сошастоп to: 


ECT, 4) = Ea + t4 + t. + (t, E t. )(1-P)/P 


Smoler also assumed that for tanks, P(K/H) = 1.0 and 


Us - O, reducing the formula to: 


ET, y? = c + a 


- t, (1-P)/P(n1m) 

гам 222405 these equations to calculate the 814 5, the 
attrition during each time step was computed using the 
Euler-Cauchy differencing equations to approximate Lan- 


Bmester's force-on-force attrition differential equations. 


Co UNIT LOCATION AND MOVEMENT 

Red and blue unit locations were handled in two dif- 
ferent ways in the original model. The blue locations 
Were left as user inputs, while the red locations had 
been pre-determined by the model builder and could not be 
altered by the user. This allowed flexibility of defensive 
positions, but because of the method of determining movement 
meuves for the attackers, this flexibility was limited. 

The method utilized for route determination was, at 
ШУ unrealistic. For each original red location, a 
straight west to east route was calculated. Each route 
was divided into 40 meter intervals, since each red unit 
was assumed to move that distance during each 10 second 
time interval. This approximated an average rate of move- 
ment ОТ 9 m.p.h. This method of route selection, although 
easy to implement, has several significant shortcomings. 


Ш (сес not allow the attacking units to utilize terrain 


d 





features during movement. Also, it permits attackers to 


Ногоо сетат that 1S, in real life, impassible» 


D. BATTLE TERMINATION 

In any combat model, adequate battle termination rules 
must be considered. Smoler utilized two criteria LO 
terminate the battle. The first of these was annihilation 
of one of the two forces. This criteria is reasonable in 
an expected value model like this and was adequately 
handled. This is not the case for the second termination 
rule. 

The second rule for termination is that the distance 
between red and blue forces becomes too small. The use 
Semis criteria for a stopping rule is fine, provided it 
is implemented properly. Іп the original model, a geographic 
center was calculated for each force. It was these center 
of mass voints that were compared to determine if the 
units were too close. This method can lead to problems, 
such as those that were encountered during some initial 
ps! runs. Іп one battle, the red units, during their 
advance, had eliminated cne entire blue unit and had suf- 
fered the loss of one of their own units. Since the red 
units are allowed to advance as much as 150 meters apart 
and blue units were set up so they were approximately 
1400 meters apart, there was considerable distance between 
units on the battlefield. Since the two units had been 


destroyed, the center of mass computation allowed the red 


ша 





units to move completely through the blue defense, while 
КЕСИ sides st111 had forces remaining to fight with. This 
reflects a problem with the center of mass computation as 


a method for determining distances between units. 


Eb. DETECTION AND FIRE ALLOCATION 
Both detection and fire allocation processes are nand- 
led well in the original model and no changes have been 
made. A detailed description of these methods is contained 
in Smoler's Thesis (Smoler,1979). However, a brief and 
general discussion of both processes will be included here, 
The detection phenomena is modeled in two ways. First, 
a non-firing detection can occur as a result of an observer's 
random search within his designated section of responsibility. 
A 30° field of vision for an observer is assumed, and the 
probability that an observer is looking in the direction 
which enables him to detect a target is computed by inte- 
р нг гс the Limicon Function over limits that are + 15° 
from the primary direction the observer is looking. The 
Паша соп Function, f(8), is the following probability den= 


ЕЕСУ function: 
£(6) = A+ Bx cos(@) 


where 
D = assigned sector width/2 


А = -B x cos(D) 


Ud 
i! 


1 / 2(sin(D)- D x cos(D)) 


© = primary direction observer is looking 
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Also, A and B are chosen such that 
D 
J sede = 1 
-D 


The second method of detection is a firing detection. 
This phenomena occurs when the following happens. If a 
mercer’ s location is within + 15° of an observer's primary 
direction of observation when he is firing, he is assumed 
to be detected and is added to the observer's target list. 
This models the detection of a firer by locating a firing 
Signature of a weapon. 

The fire allocation process is also modeled well by 
Ёо Сг, Since each- firing unit is not restricted to firing 
at one target, a fraction of each firing unit is allocated 
to each target on the firer's target list. This fraction 
is determined as a function of range to each target and of 
predetermined priorities of fire that are set so each unit's 
fire priority goes to those targets to his immediate front. 
Targets to the firer's flanks are then allocated a smaller 
percent of the available firepower. Again, the detection 
and fire allocation processes have been handled well and 


no changes to these routines will be made. 





ІІІ. MODEL CHANGES 


А. GENERAL 

As the previous section outlines, there were several 
major problems discovered in Smoler's original model. In 
order to make the model more flexible and classroom useable, 
several major changes were found to be necessary. This 
Chapter will provide a detailed description of these 
Changes, which include conditionally-deterministic attrition, 
unit locations and movement, battle termination and move» 


ment of defenders. 


“В. CONDITIONALLY-DETERMINISTIC ATTRITION 


Шимон цер Го пиггойисе a stochastic process into the 
attrition computation for this model, several options were 
available. These include the use of a Markov-process to 
determine casualties or the use of random attrition-rate 
Meeiiicients. After considering both options, and in 
keeping with the transvarent nature of this model, it was 
decided that stochastic attrition coefficients should be 
utilized, because the attrition-rate coefficient is a 
random quantity measuring a unit's fighting ability and 
can be realized before any given battle. The following 
procedure was used to develop a method for random attrition- 
mate coefficients. 


The attrition-rate coefficient, A is a measure of 


ij’ 


mie rave a firer of unit i attrits a target of uniti j. 


1. 





This also can be interpreted as a measure of the fighting 
pitty sot an elemenu or type i. it is intuitively obvious 
that this is a variable quantity that can be affected by 
many different factors, including weather, espirit-de-coros, 
previous engagements and leadership, to mention only a few. 
In a homogeneous model, one can attempt to capture these 
Mendo eifects by developing a distribution of initial 
menting unit capabilities, i.e., initial 81,5 for each 
Piet, wolnee no data is readily available that captures this 
phenomena, an attempt was made to fit a distribution to 
an intuitive feeling as to how this fighting ability varies 
Maem One Unit to the next prior to 2021 engagement. 1% 
was the author's intuitive feeling that the A, ,'8 should be 
distributed between „5 and „б with the majority of the units 
being rated between .5 and .6. 

Utilizing the above intuition, an attempt was made to 
feat a distribution to these assumed fighting levels. 


Graphically, the distribution would look like Figure 2. 


шс a UrStri bution of Initial Ajj 5 
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Initially, an attempt Was ша: тг Beta distribution 
mechs curve. The fit was fairly accurate in the middle 
range, but was not satisfactory in the tails. The same 
was true for a Normal approximation. It was finally 
determined that the best fit was a straight quadratic fit, 
where a Uniform(0,1) input variable could be used and the 
output would yield the desired random number. The quad- 


ШС fit that was used is: 


о _ 2 
Ajj = -2U 


ЕУ, 

Once a distribution for the attrition-rate coefficients 
had been derived, the implementation into the model was 
accomplished. Since it was assumed that the fighting 
КООШО ТУ Of each unit is a random quantity prior to a given 
battle, it was only necessary to obtain a realization of 
me random variable for each unit prior to the initiation 
of the battle. This realization, AY is determined by 
using a random Uniform(0,1) number and the above formula, 


A new attrition-rate coefficient for each unit is then 


computed during each time step using the equation: 


O 2 z 
А; lerne) цог 
As = 
О Tor ПЕ Ег 
where 

Ге = maximum effective range of a firer's weapon 

r = current range between firer and target 

AS IES MAIZ action ot unit's fighting ability 


12 





This equation is utilized because it varies with range, 

БИШ 15 also 15 а function of Aj yo thus creating a different 
attrition-rate curve for each unit, depending on that unit's 
fighting ability prior to the battle. Graphically, as an 
example, the attrition-rate coefficient curve for an Аз; 
of .6 and a maximum effective range, r,, of 2000 meters 


Mula look as in Figure 3. 





1000 2000 3000 


range(meters) 


I i s et ron=rate coefficient curve- for 
AS = .6 and р. = 5000 п. 

Опсе the above method for determining attrition-rate 
coefficients had been selected, it was coded and included 
as a user option in the program. Once the code had been 
implemented, several runs were made utilizing different 
random number seeds to compare battle results uSing this 
method for attrition-rate coefficient determination. These 
runs Were then compared with a run using the original 
deterministic method. The results are summarized below. 


First, the model was run using the deterministic attrit- 


ion method. For purposes of comparison, all locations, 


Ша 





movement rates, force levels and other input variables were 
med constant for all runs. Combatant attrition using the 
deterministic process as a function of time is shown in 


Figure e 


Force 
Levels 





ЕС 4O00 -500 
Time 


ши сс дева слот using deterministic attrition 


In this battle, termination occurred as a result of forces 
being too close together, with force levels at termination 
reflected in Figure L, This was then compared to the runs 
that utilized the new attrition module, The results of 
ШЕШІ ОГ these runs are in Figure 5. The IS are the 
realizations of the random variable denoting a unit's 
ИШ ла! fighting capability prior to the battle. The three 
numbers under the red and blue headings are the realizations 
for each subunit in the battle on which casuality assessments 
are computed. 

These results show, in fact, that a randomly selected 
A: does have an effect on the final outcome of the battle. 
This module now provides the user the option of selecting 
a conditionally-deterministic process that still uses 
Lanchester's equations as the basis for the 2. 


eomputation. 
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oe Force Level 


iJ 
Bus Red Result 


KEE 122127 Too close 
2205967 2601274 260 .72 Blue wins 


Figure 5. Comparison of Stochastic and Deterministic 
Battles 








С. UNIT LOCATION AND MOVEMENT 

As was pointed out in an earlier section, the unit 
locations and movement logic that Smoler used tended to 
be rigid and inflexible. The changes that have been made 
to rectify this and allow more flexibility for the user 
will be discussed here. 

Lurst. the user is now responsible for inputing all 
combatant locations, both attackers and defenders. The 
format for these inputs will be specified in Appendix A, 
the User's Manual.  Aditionally, the user is given the 
option of selecting routes for the red advance or allowing 
the model to calculate straight west to east routes as in 
phe past. To add more flexibility, the user will also be 
required to input the rate of advance (vehicle speed) for 
the attacking forces. The inputed rate of advance is 
feed in both methods of route calculation, thus controlling 


the*red attack. 





1. Straight Line Routes 


The Option to use straight line routes from west 

to east is handled in the same manner as before. The route 
is broken down into discrete distance intervals from the 
initial red location straight to the east. The only dif- 
ference is that the length of each interval is determined 
by the vehicle speed that the user inputs. The route for 
each red unit consists of 125 equal length intervals. 

2. User Determined Routes 

The option to allow the user to select routes for 

the red advance has been added to enable the user to make 
use of available terrain and to add realism to the model. 
The method for calculating these routes is straightforward 
and is discussed next. 

The user is required to input the original location 
of each red subunit and from one to ten nodes he wishes 
Beem attacking subunit to move through. This information, 
along with vehicle speed, is used to calculate route inter- 
vals that move the attacking unit through each of the 
designated nodes. A complete route would look like the one 
depicted in Figure 6. The method used to compute the routes 
is as follows. 

The straight line ground distance between the first 
two adjacent nodes, DIST, is calculated as shown in Figure 
7. The angle between the desired direction of movement 
and straight west to east movement, a, is then calculated. 


Utilizing these quantities and the distance desired to 
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move during each time step, DST, the distance to be moved 
in the x and y direction, XLN and YLN, is now computed as 
shown in Figure 7. These distances are then added to the 
coordinates of the previous interval endpoint, point С 
ne figure, to determine the coordinates of the next 
interval endpoint, point 0. This same distance is again 
added to compute the coordinates of the next endpoint, 
point E. This process is continued until the distance 
from the last endpoint computed to tne next node is less 
than DST. This process is repeated between the next two 
nodes until the entire route is completed, In order to 
insure all intervals are of equal length, the computation 
of the first interval between any two nodes must be con- 
Sidered separately, by taking into account the distance 
left over from the last computation between the previous 
two nodes. To do this, the first interval takes the rem- 
aining distance, e, and adds it to an interval length of 
DST-e for the first interval between any two nodes. This 
insures that each interval along the route is of length DST, 


Which is the required length. 


D. BATTLE TERMINATION 

As previously stated, the battle is terminated either 
Dy annihilation of one of the two forces or by forces 
Peering too Close. The annihilation criteria has not been 
Changed. However, to insure the distance criteria is 
effective, two changes have been made. 


First, the distance between units is no longer calculated 


2% 
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XLN = DST x COS(a) 


Figure 7. Route Computation 
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гоог Or mass distance. Instead, a distance is 
computed between each attacking subunit on which casualties 
are being computed that is still alive and each defending 
Eubunit that is still in the battle. If any of these 
distances becomes too small, the battle is considered to 
have moved to a "close-in" combat (hand-to-hand) that this 
program does not currently model. For this reason, the 
battle is then terminated. However, to insure the attackers 
do not pass by the flanks of remaining defenders and remain 
Outside termination distance, a check of x coordinates 

for each subunit is also made. If any attacking subunit's 
E sordinate places him past the location of the forward 
Mest blue defensive subunit still in the battle, the battle 
is also terminated. 

The criteria for being too close has been left as a 
Meer determination. it is one of the inputs that is 
described in the User's Manual, Appendix A. This allows 
ШЕШІ ТІехісііліу от breakpoint distances which also lends 
itself to the study of optimum breaking distances for 


various weapon systems on the battlefield. 


E. ALTERNATE DEFENSIVE POSITIONS 

The last addition to the model, added to increase the 
meexi bility and realism of the battle, is the option of 
Waving the defending units move to alternate positions if 
the attackers close to within a user specified range. 
The move is handled in a simple and transparent manner. 


If the user selects this cption, the defender will move 
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when the attacking forces close within the breakpoint 
Басе. When this occurs, each defensive unit that is 
still alive will move to an alternate position the user 

has selected. The duration of the move is also a user 
input. He simply specifies the number of 10 second time 
steps he wants to allow the blue forces to move. At the 
completion of the move, during which the red forces continue 
to attack, the battle will continue until one of the term- 


EN yon criteria discussed above has been met. 





IV. FUTURE ENHANCEMENT 


Although the above changes have been implemented and 
a more useable and flexible model has been created, there 
are several areas that could still be considered in future 
work. These include the introduction of heterogeneous 
forces, explicit computation of ammunition expenditure 
and artillery (indirect fire). Possible methods for 


employing these ideas are discussed below. 


A. HETEROGENEOUS FORCES 
—— The current model involves combat between homogeneous 
merces only. In other words, each opposing force is 
comprised entirely of one weapon system type. Added 
flexibility could be attained by allowing heterogeneous 
force structures on both sides. This would enable the user 
to investigate the effect of different force mixes on 
battle outcome. 

The introduction of different weapon systems within 
a single unit would require extensive changes to the att- 
n process, Although Lanchester equations for aimed- 
fire could still be utilized, casualty assesment against 
each weapon system type by each opposing weapon system 
type would have to be calculated. The total attrition of 
any particular unit would then be the sum of the damage 
assessed to each weapon system of that unit. Aditionally, 


the force level for each weapon system would be required 
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as an input as well as a separate set of hit and kill 
probabilities for each weapon system type against each 
Eve of target. 

Introduction of this feature would create a more comp- 
Jicated model, adding realism but detracting from the current 
Simple and transparent form. Since the purpose of the 
wr ent effort has been to maintain this simplicity, this 
option has not been included, but could be considered for 


шасиге work, 


В. AMMUNITION EXPENDITURE 

шде сопрас ргосезв 1s a complicated and intricate 
process. One of the hardest, and not many times attempted, 
areas to model is the logistic area. However, in the 
model currently being studied, the problem of ammunition 
expenditure could be modeled. 

Since the model is an expected value type model, it 
would not be effective to model the amount of ammunition 
expended by simply counting bullets as they are fired. 
Actually, in this type of model it would not be possible 
to count each round, because of the aggregated nature of 
each unit. The only way to model ammunition useage is to 
model it in the same way as attrition. This simply means 
that each unit would have to be given a starting level of 
ammunition on hand, a basic load, and each time the unit 
fired, the expected ammount of ammunition expended would 
be subtracted from what is on hand. The ammount expended 


would necessarily be a function of the size of the unit 
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firing, the number of targets fired at and the rate this 
Bou This can be directly correlated to the 
rones rate coerticient for the firing unit. when 
ammunition on hand reaches the zero level, the unit would 
have to be removed from the battle permanently or for some 


specified time period, to simulate resupply. 


ee ARTILLERY 

Шок armed conflier, there are more ways to inflict 
E gtres than just direct fire. liany types of indirect 
fire are utilized on the battlefield, including artillery, 
Close air support and naval gunfire. To include this type 
Ру in this model, primarily artillery, would involve 
a major effort with large program additions. 

To model artillery in an expected value model that 
already uses Lanchester equations for direct fire would 
Suggest the use of Lanchester's area-fire equations. ‘These 
equations are no more complicated to handle than tne aimed- 
fire equations that are already in use. However, other 
model considerations would have to be investigated. These 
venta include location of artillery units, whether to use 
forward observers for target location or only use pre- 
planned strikes or both and fire allocation procedures 
ШІ ре used. Also, a Significant ammount of data would need 
to be collected concerning weapon types to be used, eff- 


mn ae ranges, killing radii and accuracy data. 
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сг: 02062 20200014 ре a significant addition to 
the model. However, this too would effect the simplicity 
and transparency the model currently possesses. This 


ЕНШІСІ Пас 3159 рееп left for future consideration. 


--- u em 
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APPENDIX A 


User's Manual 


I. INTRODUCTION 


Smoler's model of land combat is a force-on-force 
combat model that utilizes Lanchester's aimed fire equations 
for casualty assessment. The battle is simulated ona 
10 x 10 KM piece of terrain representing an area east of 
Fulda, West Germany. A portion of the terrain map is 
inclosed as Figure 1. It is a computer model that is 
coded in FORTRAN and is available for use in the Y. Ro 
Church Computer Center. 

The purpose of this manual is to familiarize the user 
with the model and to provide the required Job Control 
Language (JCL) and user inputs to run the model. A 
sample listing of program output will also be provided to 


familiarize the user with the expected program output. 





Чең ителә, Тәрор 


{ 


| Sanar: 








шог GSA STRUCTURE 


Smoler's enhanced model is a computerized model, coded 
in FORTRAN, containing a main program and nine subroutines. 
To assist the user in understanding how the model operates, 
a brief description of the function of each of these major 


parts of the program will be included. 


A. MAIN PROGRAM 

The main program has several important functions. 
ШЕШ ОТ the input and output functions, except the line-of- 
sight data, are contained in the main programe Aditionally, 
the main program is used to structure all of the other 
шшисътопв during each lO second time step. Attrition, 
detection, movement and fire allocation are also handled 
in the main program. The nine subroutines provide needed 
input numbers for the above calculations. Тһе general flow 
of the program is depicted in tne flow chart shown at 


iure 2, 


ES SUBROUTINES 
There are nine subroutines included in the model. Each 
Of these subroutines performs a distinct function, eacn of 
which will be discussed below. 
2 5 oubroucine SETUP 
1015 subroutine is used to read in the terrain 
data for the Fulda Map. This terrain data will be used 


when computing line-of-sight between targets and observers, 
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INITIALIZATION: 


Read input data (terrain, initial force levels, initial 
positions, routes data, tactical parameters, etc.). 
Initialize variables and array values 


Time = Ture + 10 sec 


Update locations of attacking units 


Check line-of-siaht between any two opponent units 


NO 
LOS exists 


YES 














fire conditions are 
satisfied 


YES 
Allocate targets 
Update the accumulated detection probabilities 







Allocate fire 


Campute attrition 


Print results for current time step 


2: End of battle 


the 
YES 
Print final results 
FIGURE 2. FLOW CHART 
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“уг аз ргоупатпе a grid system for unit locations and 
movement. 
8: ошогошолие LOS 
This module, developed by Professor Jim Hartman 
Ewthe Naval Postgraduate School, computes a percent of а 
Perpeet visible to a particular observer, given the coor- 
ШЕ е5 (location) of both. This visible fraction is used 
in the detection and attrition modules in the model. 
5. Subroutine KOVER 
This subroutine is used by subroutine LOS in deter- 
mining what portion of a particular target is covered by 
the terrain between the target and observer. This number 
r used in the detection and attrition modules. 
4. Subroutine ETK 
This module computes the expected time for a given 
Тігег to kill a given target. The calculation is a function 
EE suse, time of flight for a round and hit and kill 
probabilities for the firing weapon system. It is a number 
Mat 1s used in the computation of the deterministic attrition 
Exetilcilents, А... 
Pee Subroutine STOCH 
This is used when a user has selected a stochastic 
m Lion option to compute the attrition coefficients 
tienes each time step. The calculation is a function of 
the original stochastically determined А as well as a 


mame tion of range. 





оо 12101116 11МБРА 
This subroutine is used in conjunction with tne 
r routine со compute the detection rate of a target by 
an observer given the percent of target visible to the 
Eserver, 
8. Subroutine ROUTE 
This is used to compute the route of each attacking 
red unit when the user has selected the option of inputing 
attacker routes. It calculates the coordinates of each 
interval endpoint along the route, making each interval 
length (distance moved during a 10 second time step) the 
same. The interval length is determined by the speed the 
user has selected and inputed for the current battle. 
9. Subroutine ELEV 
A subroutine that is used to calculate the terrain 
 гастоп for а given set oi X, Y coordinates. This function 
ШО е n conjunction with the LOS subroutine in comouting 


line-of-sight between observer and target. 
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е AVATLABLS OPTIONS 


The enhanced version of tnis model has been written 
to allow for maximum user flexibility while maintaining 
the simple and transparent nature of the model. To allow 
this flexibility, there have been several user options 
Mmge@r porated into the model. Each of these options, 
including user responsibilities, will be discussed here, 
(гэ Пе required input data for each being outlined in 


the next section. 


A. БКО Л ОТТО VS DETERMINISTIC ATTRITION 


02 


The user is required to specify whether he wants 
se stocnastic or deterministic attrition calculations. 
Poth methods utilize Lancnester's aimed-fire equations, 
the difference being the method of calculating the attrition- 
rate coefficients. The deterministic procedure uses the 
Bonder-Farrell method of calculating the Ai 4» while tne 
Beocnhastic method selects an initial random Ajj for each 
ШІ, апа uses this as a function of range to calculate 


the Ajj for each time steve. 


EN ATTACK ROUTES 

The second major option available to the user is 
Br zmerhod of route computation for the attacking forces. 
The user has to decide whether to allow the program to 


ШОО рше Straight west to east routes or to input the routes 





Memeesires eacid асас ле unit to follow. To select his 
own routes, a user must input the number of nodes he 
Wishes to have on each route and the coordinates of each 
of these nodes. The vrogram will tnen compute routes 
through each node. The nodes must be inputed in order 
from west to east and should not create an angle between 


Mme west to east axis and the route direction that exceeds 
O 
15°. 


Co ALTERNATE DEFENSIVE POSITIONS 

The third option the user must consider is the 
use of alternative defensive locations. This option permits 
the user, if he desires, the capability to move the defen- 
ders to alternate positions if the attackers close within 
some specified distance. This breakpoint distance is 
decided and invuted by the user and is also used as the 
Бе сапсе Гог battle termination. The alternative to moving 
the defenders is to terminate the battle when the breakpoint 


e ance is reached. 


Еэ OLHER INPUTS 

There are other inputs that are required by the 
program that the user must provide. These include force 
Sizes, weapon characteristics, unit locations and hit and 
Meee probabilities. The required formats for all inputs 


emer Ouclined in the next section. 
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IV. INPUT DATA 


All user input data is read in from cards at the start 
Шееле main program, MAIN. A brief description of each 
data element as vell as the required format for inputing 
this data is discussed next. A sample 21. has 
also been included as Figure 5, with the referenced data 
element number appended on the right. This sample data 
deck reflects the use of all data elements. However, 
not all the elements described below are required, and 
those that are needed only if a particular option is 
selected are noted with an asteric next to che data element 


number. 


Шаға Element Description and Format 





Ж О Ше л О атап Option and random 
number seed. A lis entered for deter- 
Шин стоте на О or stochastic attrition. 
This number is followed by a five digit 
remem Wumoer сес roma for this card 
werd, 

2. Number of blue units (NBU) and number 
of red units (NRU). Both of these 
elements are two digit integers. Card 
о Ее ра ле пе. 

Ё» Effective Weapon ranges listed in the 


following order. Red minimum, red 





6. 


p." 


maximum, blue minimum and blue maximum. 
СОТ Ко тзг їз L(F6,1,1X). 
Red force levels for each unit. 
ar Ола as 85.6, LX. The card will 
contain the number of entries that equals 
Blew nunber em redwunits (WRU) in data 
element 2. 
Tooo Toute Computation desired followed 
by vehicle speed. Use a 1 for user 
determined routes or a О Рог progran 
determined routes. A one digit entry will 
designate desired speed as follows: 

1 for 9 mph 

OLE mph 

bur lo npa 

ог То mon 
Card format for these elements is le, 
мо Те. 
E G Oo r Les 10 each initial геа 
location, One card is needed for each 
amu tavern Tormat Poel, 1X, F6.1. 
If the user has selected to enter his 
own routes in data element 5, the 
following is required for each route. 

а. The number of nodes in the route 
(frommol to 10). This card's format 


n SUO. 


ШІ 





9. 


Uo. 


ШІ. 


be X, Y coordinates of each node 

along the route in order from west 

Bow Опе сава Рог each node in 

що иа Вой, Пи 56,1. 
тӘ асе пог ol ue location, force 
Ше) топ васп unit. ока пагу search direction 
for that unit and desired search width. 
Force level must be between 1.0 and 3.0, 
search direction between 135° and ӘНЕ Б 
and search width from EOS to 1209. One 
гс Ол Бе ИЕ with format F6.1, 
E 990120715, 5, 1Х, 13. 
DECIS Me is to move to alternate 
defensive positions, breakpoint distance, 
and number of lO second time intervals 
allowed for the move. A I for no move 
ТІН сес! tO move, Kormat for thas 
ОБО . moll IX. Ie. 
Alvernate blue position X and Y coordinates 
if move is desired (data elenent 9). One 
ud ropSccche location with format F6.L, 
т.т 
Red weapon system hit and kill probab- 
ilities. The probabilities are entered 
шог атг 01-20 опе сага for each 
range band. Four probabilities for each 


range are probability of a 75 Found nat, 


ТІП 





probability of а hit given a previous round 
o OD ability Or a піс given a previous 
FoRo miss, and probability of a kill 
given a hit. The 6 range bands are from 
ЖЕСТ 50077571909, 1000 to 1500, 
ШОО О ШЕ 2000322000 te 2500 and 2500 To 5000. 
Each of these six cards has the format 
DX oe. 
ша. Same as item ll, except the vrobabilities 
for the blue weapon systems are entered. 
ос 0020-0121 format h(Fh.2, IX). 
вне formats listed above, a lX is a space, an IX is an 
integer of x digits, and a Fx.y is a real number of length 
И У digits to the right of the decimal. For example, 
НЕ СОПІП be a 25, an T6.l could be 3487.4. A format 
Containing these specifications could be 12, 1X, F6.1. 
The data card would then look like 25 5487.4, witn the 25 


wenning in card column onee 


дг 





card 1 2 5 L, fent 
column 1234567890123456789012345678901234567890 number 
1 28945 1 
05205 с 
ОО 100040 3 
ЖОО ер L 
1 2 5 
7500. 07750020 
2400.0 5900.0 6 
A OOOO 
02 


4900.0 7200.0 

7100.0 6700.0 

бг 

БЕОСЬС 500050 7 
7000.0 6500.0 

on 

8900.07 550050 

27,002077200207 5201907120 

2 21201-:5002022,02120-120 9 


7700075290207 7.071°807 1720 
A A tw 


22705005071 9 
7000.0 7200.0 
7100.0 6500.0 10 


7100.0 5800.0 


On 202970875 
0.80 0.80 0.75 
222 275 10.70 
О ООО Сб О6О 
0.45 0.50 0.50 
ео 0.20 0.20 


Fi re 5, 


0.70 
0.70 
0.65 
0.55 
0.22 
0.20 


sample Data Deck 


j] 





ward 


column 


1 2 - 
1234567890123456789012345678901234567890 


0.60 


0.85 
0.80 
0.75 
0.60 
0.40 


0.70 
0.90 
0.85 
0.80 
0.70 
0.145 


0.65 
0.85 
0.85 
0.75 
0.65 
0.40 


ү илге 5, 


222 
0.90 
0.80 
0.70 
0.65 
0.50 


Sample Data Deck (cont) 


Ll 


даса 
element 
number 


ШЕ! 





У. REQUIRED JCL 


The model and line-of-sight data are residing on 
permanent disk in the W. R. Church Computer Center. Іп 
order to exercise tne model, the required job control 
language (JCL) is illustrated in Figure 4. Due to the 
CPU time required for execution, a time parameter of 


TIME = 2 must be used on the job card. 


//JOB CARD , TIME=2 

// EXEC FORTCIG 

//FORT.SYSIN DD UNIT=2314,VOL=SER=PATOO2,DISP-SHR, 

// DSN=S1360.SMOLER. PGM, DC B=(RECFM=FB, IREC L=80, BLKSIZE= 3200) 
2» 

//GO.FTO8FOO1 DD UNIT=2314,VOL=SER=PATOO2,DISP=SHR, 

77 DSN=S1360.LOSDATA,DCB=(RECFH=FB, LRECI=80,B1KSIZE=3200) 

/» 

//GO.SYSIN DD * 


USER DATA DECK 


/* 


Figure L, Required JCL 
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Е БЕСТЕ OUP PUT 


Once the model has been executed, the user can expect 
a following output: 
i -rogram listing of the model, 
e- A summary of the initial battle conditions, 
ОСОРИ аа пе ишосастопе and options selected. 
3 - A summary of battle conditions after each 10 
second time step, including unit locations, force 
levels, unit status, percentage lost and targets 
on each unit's target list. The status will 
show one of the following: 
О = Unit alive, not firing 
l - Unit alive and firing 
2 - Unit killed 
3 - Unit moving 
1 = A statement of the reason for battle termination. 
A sample of the initial battle summary is at Figure 5 and 


the time step summary is at Figure 6. 
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INITIAL BATTLE INFORMATION 


UNIT X Y FORCE LEVEL 
1 1900.0 7800.0 860 
2 ШЕШЕСІ 6500,0 S 
3 2020 1000209 6.0 
L 25 01:02 ООО 2.0 
D "oo o SIC 5G 250 
6 6200.0 5150.0 p 


ATTE TONT TS STOCHASTIC 
ROUTES DETERMINED BY USER 

poe nO tomes 15 15.0 
BREAKPOINT DISTANCE IS 1500.0 


BLUE WILL MOVE TO ALTERNATE POSITIONS 
ALTERNATE POSITIONS ARE: 


UNIT 
L 


2 
6 


X 
9500.0 
9800.0 
9200.0 


ү 
770030 
6200.0 
5800.0 


RED KILL PROBABILITIES 


RANGE 
500 
1000 
1500 
2000 
2900 
5000 


| 
0.95 
0.80 
0.65 
0.40 
0.30 
0.20 


BLUE KILL 


RANGE 
500 
1000 
1500 
2000 
2500 
5000 


Figure 


P 
0.80 
0.80 
O. 79 
0.65 
0.45 
0.42 


De 


PHH 
0.92 
0.85 
0.70 
0.50 
0.40 
0.30 


PHM 
0.95 
221! 
0.70 
0.22 
0.40 
5/22 


РКН 
1.00 
0.95 
0.90 
0.75 
0.75 
0.60 


PROBABILITIES 


PHH 
0.80 
0.80 
0.7/0 
OG 
0.22 
0.22 
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РНМ 
0.80 
0.80 
0.70 
0.65 
0.45 
0.45 


PKH 
1.00 
0.95 
0.90 
0.85 
0.70 
Onc © 


Initial Battle Summary 





TIME = 480 SECONDS 


ENLT 


ONV FE XA Го E 


Х 
4255.2 
4852.53 
4502.6 
5800.0 
7400.0 
6200.0 


Y 
229952 
2799.9 
4279.0 
7700.0 
6300.0 
2150.0 


паке о 


FORCE LEVEL 


N— wo O O O 


Time Step Summary 
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STATUS  LOST-PCT 


TARGETS 





о соз Ф лас шом © 


г 


Сэ e» o o 


о г 


сэ 


APPENDIX 3 


Program Listing 


COMMON /GRP1/ ІРАОІА (6) , І15ЕСМО (6) , MVTDIR (6) , X (6) , T (6) , SPD (G1 
COMMON /GRP2/ TR(2)1,T1 (2), TH(2) , TM (2) , TF1 (2), TF2 (2) , TF3 (2) , 

1P (2,6) , PHH (2,60 , PHA (2, 89) , PKH (2,6) , TF (2) 

COMMON /GRP3/ NBU,NRU,FL (6) , FO (6) , NO1 (3 , XIC(3,200), T1C (3,2001 , 
1101R (3,200) ,AVSP, 1SPD 

1, IUSTRT (6) , 11 (6) , LOST (6,6) ,VISFRA, VISFAB,SIZETK, 
ISIZETW,NT (6) , NF (60) , SRF, DISMRX, 

INLOSC (6,6) , VISFR (6, 6) , RHMINTK, RMXTK, RMINTH, RMXTH, OP, TOHFR, TNKFR, 
1PTT (3, 3) , RF, P0R (6,6) , RPOR (6,6) , LOR (6, 6B) , NR (O) , OFL (6) , POL (61 
COMMON /GRPu/ TPOL (6) , 8.00 (6,6) , 0 (6, 6) 

COMMON /GRPS/ LOT (6,6) , ROT (6,6) 

COMMON /HILLS/ XC (100) , YC (100) , PEAK (100) ,SX (100) ,SY (100) , RHO (100) 
COMMON /HILLS/ SCALE (100) , TNORHO (100) , TWOSCL (100) „BASE 

COMMON /HILLS/ NHILLS 

COMMON /COVER/ CXC (150) , CTC (150) , CPERK (190) , CPXX (190) , CPYY (150) 
COMMON /COYER/ CPXY (150) , NCVELS 

COMMON /COUNTR/ KH,KHMH,KV,KN,KGRS, KELL, KINT 

COMMON /GRID/ LST (10,10) , NHL (10, 10) ,LISTH (450) , KHREP (100) , KTREP 
COMMON /GRIO/ LSTC(10,10) , NC (10,10) , LISTC (400) , KCREP (130) 
COMMON /GRP6/ ALPHA (6) 

COMMON /GRP7/ XA (6) ,YA (6) , 1MOVE (6) 


INITIALIZATICN. 


CALL OVFLON 
В: 0.0 

А «0.0 

МР«0 
PAI=22.0/7.0 
ZL».00001 
CRLL SETUP 


AERD TERRAIN DATA FOR LINE OF SIGHT 
CHECK FOR STOCHASTIC OR DETERMINISTIC ATTRITION 


READ(S,130) ITRIT,1S 
130 FORMAT (11,1X,15) 

DC 132 1«1,6 

CALL RANDOM (15, TRAN, 1) 

RLPHA (1) » (2-2. *WYTRRN*x2) * (2. *XTRRN*.3) 
132 CONTINUE 
READ IN NUMBER OF BLUE AND RED UNITS 


ВЕВО (5,200) NBU,NRU 
200 FORMAT (12,1X,12) 


INITIALIZE WEAPON SIZES 


512ЕТК-2.5 
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51 
52 
33 
su 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 


72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
au 


85 
87 
88 
88 
90 
91 


93 
94 
95 


97 
98 
99 
100 


om 


о 


сэ 


гэ 


SIZETN=2.5 


READ IN EFFECTIVE WEAPON RANGES 


102 


READ (5,102) RMINTK, RMXTK, RMINTW, RMXTW 
FORMAT (F6.1,1X,F6.1,1X,F6.1,1X,F6.1,1X) 


INITIALIZE PM,RF,TOWFR, TNKFR AND NOD 


101 


111 


PM=.352 

RF=.5 
TOWFR=.03 
TNKF R=. 1 
NOD=2 

06 101 ]=1,NRU 
NOI (1) 2125 
CONTINUE 
KzNRU*1 

Ls NRU*NBU 

DO 111 l=1,L*“ 
11 (1) =0 
CONTINUE 


READ IN FORCE LEVELS OF EACH RED UNIT 


103 


READ (5,103) (РЁ (1), 1=1, МАО) 
FORMAT (3 (F3.1,1X)) 


CHECK FOR TYPE OF ROUTE DETERMINITION 


106 


READ (5,106) IRTE,ISPD 
FORMAT (11,1X,11) 

12 1 5Р0.Е0.1} ЯУЗР=9.0 
IF (1SPD.E0.1) DST=40.232 
IF (ISPD.E0.2) AVSP=12.0 
IF(1SPD.EQ0.2) DST=53.643 
IF (1SPD.E0.3) AVSP=15,0 
ІҒ(15Р0.Е0.3) 05Т-67.053 
ІҒ (15Р0.Е0.Ч4) AVSP=18.0 
1Е (15Р0.Е0.Ч) DST=80.463 


READ IN INITIAL RED LOCATIONS 


6 


DO 6 1=1,NRU 
ПЕН 107) Х1о11,1),116(1,1) 


107 FORMAT (F6.1,1X,F6.1) 


CONTINUE 

ЇР (1НТЕ.ЕО.1) 60 10 108 
DO 2 1=1,NRU 

DO 2 J=2,125 


p 





101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
115 
11Ң 
115 
115 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
123 
130 
131 
132 
133 
13Ч 
135 
136 
137 
138 
139 
140 
141 
142 
183 
144 
145 
146 
147 
148 
149 
150 


108 
109 


Mete be (11, J=1) +DSTx (9-1) 
X1C (1,3) =X1C (1, J-1) +DSTx (3-1) 
IDIR(1,J) =0 

CONTINUE 

60 10 109 

CALL ROUTE 

SUMRO=0.0 

00 3 121, МАЦ 

РО (1) FLI) 

SUMRO=SUMRO+FO (I) 

X tI) =X1C (1,1) 

Y (1) =YIC (1,1) 
MVTDIR (1) =IDIR (1,1) 

SPD (1) »AVSP 

IUSTRT t1) 30 
IPRDIR (1) =IDIR (1,1) 
ISECHD 11) =120 

МЕ (1) #1 

11 (1) =] 

CONTINUE 


C READ IN INITIAL BLUE LOCATIONS 


3 
C 
C 
104 
y 
С 


SUMBO=0.0 

DO u IsK,L 

READ (5,109) X(D,Y(D),FL (1), IPRDIR (1) , 1SECHO (1) 
NOMAS FS.1, 1, F3.1,1X, 13,1%, 13) 
FO (1) =FL (1) 

SUMBO=SUMBO+FO (1) 

MVTDIR (1) «0 

SPD (1) =0.0 

IUSTAT (1) =0 

IMOVE (1) =0 

CONTINUE 


С CHECK FOR ALTERNATE BLUE POSITIONS AND READ IN IF WANTED 


400 


402 
401 


READ (5,400) IALT,BREAK, ITEM 
FORMAT (11,1X,F6.1,1X,12) 
IF (IALT.E0.1) GO TO 401 
DO 402 IsK,L 

АЕЯО (5, 107) ХА (1) „ТЯ (1) 
CONTINUE 

DELT=10. 

ТЯ (1) =20. 

T1 (1) =8, 

TH (1) x8. 

TM (1) 210. 

TF1t1) 21. 

TF2(1)21. 


pu 





151 
152 
153 
154 
155 
156 
152 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 


E 


С 


ТЕЗ), 
TA (2) «20. 
ТІ (2) =8. 
TH (2) =8. 
ТМ (2) «15. 
ТЕ1 (2) «10. 
TF2 (2) »12. 
ТЕЗ (21 =15. 


READ IN HIT AND KILL PROBABILITIES 


515 
S14 


31 


00 5 141,2 

DO 514 1-1,6 
RERD (5,515) P (I,jJ) ,PHH(1,JjJ) ,PRPM(I,J) ,PKH(1,4) 
FORMAT (4 (FY.2,1X)) 
CONTINUE 
CONTINUE 

PTT (1,1) «1.0 
PTT (1,2) 30.8 
PTT (2,2) =0.2 
PTT (1,3) 30.8 
PTT (2,3) 30.15 
РТТ (3,3) =0. 05 
00 31 1=1, МАО 
00 31 J=K,L 
NLOSC (I, J) =0 
NLOSC (J, 1) 20 

0 (1,3) =1.0 
Q(J,1) 21.0 
VISFR(I, J) =0.0 
VISFR(J,1) 20.0 
CONTINUE 

IC=1 


C PRINT INITIAL BATTLE INFORMATION 


С 


599 


600 


602 


601 


604 


603 
606 


WRITE (6,599) 

FORMAT (*1*,1X, "INITIAL BATTLE INFORMATION”) 
WRITE (6,600) 

FORMRT (/1X, 'UNIT', 7X, 'X',8X, "Y ',uX, 'FORCE LEVEL *) 
00 601 Is1,L 

НАТТЕ (6,602) 1,X(1),Y(1),Ft (1) 

FURAN (IX, 13,35X%,F7.1,2X,F 7.1, 7X.F 3.1) 
CONTINUE 

IF (1TRIT.EQ.1) 60 TO 603 

WRITE (6,604) 

FORMAT (/1X, “ATTRITION 15 STOCHASTIC °/) 

GO TO 605 

WRITE (6,606) 

FORMAT (/1X, ‘ATTRITION IS OETERMINISTIC*/) 





201 
202 
203 
204 
205 
206 
207 
208 
203 
210 
211 
212 
213 
214 
215 
216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
езг 
233 
234 
235 
236 
237 
238 
239 
240 
241 
242 
243 
244 
245 
245 
247 
248 
249 
250 


сэ 


605 IF(IRTE.EQ.O) GO TO 607 
WRITE (6,608) 
608 FORMRT (1X, 'ROUTES DETERMINED BY USER'/) 
607 WRITE (6,609) AVSP 
609 FORMAT (1X, 'RED VEHICLE SPEED IS ',Fu.1/) 
WRITE (6,610) BREAK 
610 FORMAT (1X, ‘BREAKPOINT DISTANCE IS “,Ғ6.1/) 
IF (IRLT.EQ.O) GO TO 615 
WRITE (6,620) 
620 FORMAT (1X, ‘BLUE WILL NOT MOVE TO ALTERNATE POSITIONS’/) 
60 78 625 
615 WRITE (6,630) 
630 FORMAT (1X, “BLUE WILL MOVE TO ALTERNATE POSITIONS ‘/1X, 
1 'ALTERNATE POSITIONS ARE: */1X, *UNIT*,5X, *X",8X,'Y”) 
00 635 I=K,L 
WRITE (6,640) 1,XA(1) YA (1) 
640 FORMAT (1X,13,3X,F7.1,2X,F7.1) 
635 CONTINUE 
625 IRAN=500 
WRITE (6, 645) 
645 FORMAT (/uX, 'RED KILL PROBRBILITIES'/1X, 'RRNGE',UX, 'P', 
1uX, 'PHH', 3X, 'PHM ',3X, `PKH`) 
DO 650 121,6 
NRITE (6,655) IRAN,P (1,1) ,PHH (1,1) ,PHM (1,1), ,PKH (1,1) 
6S5 FORMAT (2X, I4, 4 (2X, F4.2)) 
I RAN= I RAN+500 
650 CONTINUE 
IRAN=500 
WRITE (6,660) 
660 FORMAT (/¥X, “BLUE KILL PROBABILITIES ‘*/1X, ‘RANGE *,UX, P °, 
lie PHH Y 3SX, PHH, 3X, `PKH`) 
06 665 121,6 
НАІТЕ (6, 655) IRAN,P (2,1) ,PHH (2,1) ,PHM (2,1) ,PKH (2, 1) 
IRAN=IRAN+SOO 
665 CONTINUE 
WRITE (6,670) 
670 FORMRT('1',10X, 'BRTTLE BEGINS'//) 


UPDATE LOCATION OF RED UNITS. 


DISMAX=5000.0 
67 06 9 1-1, НАЦ 
IF (JUSTAT (1) .E0.2) GOTO 9 
ЇР (105181 (1) .ЕО.0) 6010 76 
NF (I) 3NF (1) *1 
IF (NF (1).LT.NÓO) GOTO 9 
NF (1) 21 
75 00 11 J = 1, NRU 
IF (J .EQ. 1) 60 TO 11 
IF (IUSTAT(J) .EQ. 2) GO 10 11 





251 
252 
253 
258 
255 
256 
257 
258 
259 
260 
261 
262 
263 
26Ч 
265 
266 
267 
268 
263 
270 
271 
272 
273 
274 
275 
276 
277 
278 
273 
280 
281 
282 
283 
284 
285 
286 
287 
288 
289 
290 
291 
292 
293 
294 
295 
296 
297 
298 
299 
300 


сэ 


11 


17 


133 


18 


DIST a X1(1) - X(J) 

IF (DIST .GT. DISMAX) GO TO 9 
CONTINUE 

1111) = 11 (1) + 1 

K?7211 t1) 

X (1) «xX1C (1,K7) 

Y (1) zYIC (1,K7) 
MVTDIR 11) =IDIR 11,K7) 
ЇРВО1Н (1) 4101Н(1,К7) 

CONTINUE 


LINE--GF-SIGHT CHECK BETWEEN UNITS AND TARGETS SELECTION 


DO 17 JaK,L 

NT (J) 20 

CONTINUE 

00 12 151, МАУ 

NT (1) =0 

IF (IUSTAT (1) .£0.2) GOTO 12 

DO 16 JsK,L 

IF (IUSTRT(J).EQ. 2.0R. IUSTRT(J).EQ.3) GO TO 168 
XX1sX (1) 

YY12Y t1) 

CALL ELEVI1XX1,YY1,TMACI) 

XX2sX (J) 

YYƏsT (J) 

CALL ELEV(XX2,YY2, TMACJ) 

LATOB=1 

LBTOA=1 

CALL LOS (ХХ1,ҮТІ,ТМАСІ,0.0,512ЕТк, ХХ2, ҮҮ2, ТМАС.),0.0,512ЕТИ, 
ILATOB,LBTOA, VISFRA, VISFRB) 

VISFR(1,J) =VISFRA 

VISFR(J,1) =VISFRB 

IF (VISFRA.GT.ZL) GOTO 18 

LOST (1, J) «0 

LOST (J, 1) =Q 

NLOSC (1, J) =NLOSC (1, J) +1 

NLOSC (J, 1) =NLOSC (1, J) 

GOTO 16 

LOST (1, J) 1 

LOST (J, 1) »1 

NLOSC (1, J) =O 

NLOSC (J, 1) 30 

RANGE=SQRT ( (X (1) -X (J) ) <=2+ (Y (1) - Y (3)) ««2) 
IF (RANGE.LT.RMINTK.OR.RANGE.GT.RMXTK) GOTO 20 
IF(Q(I, J).EQ.1.0) GOTO 20 

IUSTAT (1) =1 

NT (1) sNT (1) *1 

M=NT (1) 

LOT (1,M) =J 





301 
302 
303 
304 
305 
306 
30? 
308 
309 
310 
311 
312 
313 
314 
315 
316 
317 
318 
319 
320 
321 
322 
323 
324 
325 
326 
327 
328 
329 
330 
331 
Jac 
333 
334 
335 
336 
337 
338 
339 
340 
341 
342 
343 
344 
345 
346 
347 
348 
349 
350 


€» 


20 


16 
le 


25 


79 


37 


23 


ROT (1,4) =RANGE 

ТЕ (м. ЕШ.1) 5010 20 

CALL SORT (1,M) 

IF (RRNGE. LT. RMINTH.OR.RRNGE.GT.RMXTH) GOTO 16 
IF (Q(J1,1).€0.1.0) 6078 16 
IUSTAT (J) #1 

NT (J) NT CJ) *1 

Ma NT (J) 

LOT (J, M) 31 

ROT (J, M) =RANGE 

IF (M.EQ.1) GOTO 16 

CALL SORT (J, M) 

CONTINUE 

CONTINUE 

00 25 151, НАЦ 

IF CJUSTAT (1) .EQ.2) GOTO 25 
ЇР (МТ (1) .МЕ.0) 6010 25 
IUSTRT (1) =0 

NF (1) =0 

CONT INUE 

00 78 JsK,L 

IF CIUSTAT (J) .EQ.2.0R. IUSTAT (J) .EQ.3) GO TO 79 
IF (NT (J) .£Q.0) IUSTAT (J) =0 
CONTINUE 


UPDATE OF THE ACCUMULATED DETECTION PROBABILITIES. 


18841 

188=NRU 

ICC=K 

10D=L 

FR«TONFR 

OP РМ 

09 14 14188, 188 

1F (JUSTAT (1) .E£0.2.0R. IUSTAT (1) .E0.3) GO TO IM 
DO 19 J-1CC, 100 

PROP=0.0 

IF (IUSTRAT (J) . EQ. 2. OR. IUSTAT(J).EQ. 3) GO TO 19 
6100 (1, 70 20 (1, .)) 

IF (LOST (1,J3).£0.0) GOTO 15 

1F (NT (1) .6T.0) GOTO 22 

PCTVIS=VISFRA (4,1) 

CALL LAMDA (1, J,PCTVIS, DETRAT,PSUBK) 
QV=EXP (- (FL (1) x«DETRRT«OPXDELTXFL (J))) 
IF (NT (J).GT.0) GOTO 23 

Q(1, J) =Q (1, J) xQV 

GOTO 19 

QP= (1.0-PSUBK) =» (FR»DEL TFL (1) FL (3) ) 
Q11,J)=0 (1, J) » (Q0V+QP-QV»0QP) 

GOTO 19 


29 





351 
852 
353 
354 
355 
356 
357 
358 
859 
360 
361 
362 
363 
364 
365 
366 
367 
368 
369 
370 
87! 
372 
473 
374 
375 
376 
82? 
378 
379 
380 
381 
382 
383 
384 
385 
386 
387 
388 
389 
390 
391 
392 
393 
394 
395 
396 
397 
398 
399 
400 


С 


22 


32 
35 


77 
33 


24 


36 


34 


19 
14 


38 
28 


NS=NT (7) 

DO 24 11=1,NS 

K1sLOT (1,11) 
ANG1=ATAN2 (Y (K1) -Y (1) ,X (K1) -X (1)) 
ANG2=ATAN2 (Y (3) -Y (1) ,X (4) 2X QD) 
IF ((ANG1»ANG2) .GE.0.0) 6010 77 
IF (ANG2.LT.0.0) GOTO 32 
ANG=2xPAI+ANG1-ANG2 

GOTO 35 

ANG=2*PAI+ANG2-ANGI 

IF (ANG.GT.PAI) ANG=2xPAI-ANG 
GOTO 33 

ANG=ABS (ANG2-ANGI) 
AA=15.0xPA1/180.0 

IF (ANG.GT.AA) GOTO 2y 
РАОР=РАОР*+РТТ (11, №5} 

CONTINUE 

IF (PROP.EQ.0.0) GOTO 34 

IF (NT (J) .GT.0) GOTO 36 

CALL LAMDA(I,J,PCTVIS,DETRAT, PSUBK) 
DE TRAT=DETRAT RF 

QV=EXP (- (FL (1) xPROPxOETRAT=DELTFL (J) )) 
О (1,2) =0 (1, J) xQV 

GOTO 19 

Q(1,J) =0.0 

GOTO 19 

IF (1AA.EQ.1) GOTO 19 

Qt1,J) 21.0 

GOTO 19 

IF (NLOSC (1,J).LE.3) GOTO 19 
Q(1,J) 21.0 

CONTINUE 

CONTINUE 

1F (IAA.EQ.K) GOTO 38 

FR=TNKFR 

JAA=K 

IBB=L 

ІССг1 

10DzNRU 

0Р=1.0 

GOTO 37 


FIRE ALLOCATION. 


DO 28 I=1,L 

МА (1) 20 

DO 26 I=1,L 

ЇР (105ТАТ(1).Е0.2.08. 105ТАТ (1) .Е0.3! 060 Т0 26 
ЇР (МТ(1) .Ейд.0) 6010 26 

DO 27 J=1,3 


56 





401 
402 
403 
404 
405 
406 
407 
408 
409 
410 
411 
412 
413 
414 
415 
416 
417 
418 
419 
420 
421 
422 
423 
424 
425 
426 
427 
428 
429 
430 
431 
432 
433 
434 
435 
436 
437 
438 
439 
440 
441 
que 
443 
uuu 
uus 
446 
uu? 
qu8 
uug 
450 


el 


30 


uu 


29 


78 


АРОЯ (І, Ј) =0. 0 

CONTINUE 

IF (NT (1) .£0.1) GOTO 78 

IF (NT (1) .E0.2) GOTO 29 

NOT=3 

MM1=L OT (1,1) 

MM2=LOT (1,2) 

MM3=LOT (1, 3) 

PROB= (1.0-0 (1, MM1)) xQ (1,MM2) <Q (1,MM3) 
APOA (1,1) =APOA (1,1) +PTT (1,1) »=PROB 
PROBsQ(I,MMI)*(1.0-Q(1, MM2)) xQ(1, MM3) 
ЯРОВ (1,2) -ЯРОВ (1,2) 9РТТ (1, 1) «РВОВ 
PROB=0 (1,MM1) xQ (1,MM2) « (1.0-Q (1,MM3)) 
APOA (1,3) =APOA (1,3) +PTT (1,1) «PROB 
PROB= (1.0-0 (1,MM1)) « (1.0-Q (1,MM2)) «0 (1, ММЗ) 
APOA (1,1) =APOA (1, 1) +PTT (1, 2) «PROB 
APOA (1, 2) =APÓA (1,2) +PTT (2,2) *PROB 
PROB= (1.0-Q (1,MM1)) 0 (1, MM2) x (1.0-Q (1,MM3) ) 
АРОЯ (1, 1) =APOA (1, 1) +PTT (1,2) -PROB 
ВРӨОВ (1,3) =APOA (1,3) +PTT (2, 2) «PROB 
PROB=0 (1,MM1) x (1.0-Q(1,MM2)) < (1.0-Q (1,MM3) ) 
APOA (1,2) =APOA (1,2) +PTT (1,2) »PROB 
APOA (1,3) =APOA (1,3) +PTT (2, 2) «РАСОВ 
PROBs (1.0-Qt1,M4M1))«(1.0-Q(1,MM2))*(1.0-Q(CI, MM3)) 
ЯРОЯ (1,1) =ЯРОЯ (1,1) +РТТ (1, 3) «РВОВ 
APOR (1,2) =APOR (1,2) +PTT (2,3) «PROB 
АРОА (1, 3) =АРОЯ (1, 3) +РТТ (3, 3) «РАВ 

DO ЧЧ J=1,NOT 

KK=LOT (1, J) 

МА (КК) МА (КК) + 1 

IN=NA (KK) 

LOA (KK, IN) =] 

POA (KK, IN) =RPOA (1, J) 

CONTINUE 

GOTO 26 

NOT=2 

MM1=LOT (1,1) 

MM2=L OT (1,2) 

РВОВ- (1.0-0(1,ММ1)) х0 (1, ММ2) 

АРОЯ (1, 1) =АРОА (1, 1) +РТТ (1, 1) «РАОВ 
PROB=0 (1, MM1) x (1.0-0 (1,MHM2) ) 

APOA (1,2) =APOA (1,2) +PTT (1,1) ХРАОВ 
PROBs (1.0-Q (I, MM1))* (1.0-Q (1,MM2)) 
APOA (1, 1) =APOA (1, 1) +PTT (1,2) <PROB 
APOA (1,2) =APOA (1,2) РТТ (2, 2) хРВОВ 
6070 30 

NOT=1 

MM1=LOT (1,1) 

PROB=1.0-Q(],4M1) 

APOA (1,1) =APOA (1,1)+PTT (1,1) <PROB 


240 





ys] 
ц52 
453 
usu 
455 
456 
ys? 
458 
459 
460 
461 
462 
463 
464 
465 
466 
467 
468 
469 
470 
y71 
472 
473 
474 
475 
476 
477 
478 
479 
480 
48] 
482 
483 
чау 
485 
486 
487 
488 
483 
490 
491 
492 
493 
чоч 
495 
496 
497 
498 
499 
500 


сэ 


e 


26 


Че 
43 


131 


5000 
41 
4? 


46 


60 


40 


58 


6010 30 
СӨМТІМШЕ 


ATTRITION. 


SUMR=0.0 

SUMB=0.0 

DO 40 1=1,L 

IF (IUSTAT (1) .£0.2.0R. IUSTAT (1).£0.3) CO TO 40 

M6=NA (1) 

SUM=0.0 

ЇР (Мб.Ей.0) 6010 Ч7 

06 ЧІ Ј=1, М6 

M7=LOA (1, J) 

IF (M7.LT.K) GOTO 42 

ITYPE=e 

GOTO 43 

ITYPE=1 

RANGE=SQRT ( (X (1) -X (M7) ) «x2* (Y (1) - Y (M7) ) же) 

IF (ITRIT.£Q.1) GO TO 131 

CALL STOCH(ITYPE,RANGE,AJI) 

60 TO 5000 

CALL ETK(ITYPE,RANGE,T) 

AJ1=1.0/T 
SUMsSUM*RJIXFL (M7) x*POR (I, J) xDELT 

CONTINUE 

OFL (1) =FL (1) 

FL (1) =FL (1) -SUM 

1F (FL (1) .GT.ZL) GOTO 46 

FL (1) 20.0 

IUSTAT (1) =2 

IF (I.LT.K) GOTO 60 

SUMB=SUMB+FL (1) 

TPOL (1) = (FO (1) -FL (1) ) УРО (1) 

GO TO 40 

SUMR=SUMR+FL (1) 

TPOL (1) = (FO (1) -FL (1)) /F0 (1) 

CONTINUE 


PRINT AND CHECK FOR BATTLE DETERMINATION. 


ITIME=1C>»10 

DO 5? IsK,L 

IF(IUSTAT(1I).EQ.2) GO TO 57 

DO 58 J=1,NRU 

IF (IUSTAT (.)) .Е0.2) 60 170 58 

CHE CK=xX (1) -X (J) 

AVD=SQRT ( (X (1) -X (3) ) »«2+ (Y (1) -¥ (5) ) 2) 

IF (AVD.LT.BREAK.OR.CHECK.LT.S0.) CO TO 250 
CONTINUE 


58 





501 
502 
503 
504 
505 
506 
507 
508 
509 
510 
511 
512 
913 
514 
515 
516 
517 
518 
519 
520 
521 
522 
523 
524 
525 
526 
527 
528 
529 
330 
591 
532 
533 
534 
535 
535 
53? 
538 
539 
540 
541 
542 
543 
Syy 
S45 
546 
54? 
548 
549 
550 


С 


сэ 


57 CONTINUE 
GO 10 99 


COMPLETE BLUE MOVE 


250 DO 251 1=K,L 
]1F (IRLT. EQ. 1.0R. IMOVE(ID).EQ. ITEM) GO TO 6000 
IF (IUSTAT (1) .EQ.0) IUSTRT (1) 23 
IMOVE (1) =IMOVE (1) +1 
ТЕ СІНШҮЕГІЛ ЕТ. ІТЕМ) 00 TO 251 
X (1) =XA (1) 
Y (1) =YA (1) 
IF C1USTAT (1) .EQ.3) JUSTAT (1) =0 
251 CONTINUE 
99 WRITE (6,112) ITIME 
112 FORMAT (////1X, ‘TIME = ",14,1Х, “SECONDS ‘//) 
WRITE (6, 113) 
113 FORMAT (1X, “UNIT*,5X, °X*, 8X, °¥*°,5X, “FORCE LEVEL*,2X, "STATUS', 
12X, 'LOST-PCT ', 2X, “TARGETS *) 
DO 59 I=1,L 
N6=NT (1) 
IF (N6.NE.0) GO TO 48 
WRITE (6,264) 1,X (1) ,Y (1) ,FL C1), TUSTAT CI), TPOL (1) 
264 FORMAT (3X,11,3X,F7.1,2X,F7.1,6X,F3.1,9X,11,6X,F5.3) 
50 TO 59 
48 НА1ТЕ (6,114) I,X(IJ, Y G) ,FLIGD,IUSTRTID , TPOL (OD, 
1 (LOT (1, J) , J21, N6) 
НИВЕЛИР 7 12X,F7.1,0X,P3.1,23X,11,6X,F5.3,3X,3(11,1X)] 
59 CONTINUE 


CHECK FOR BATTLE DETERMINATION. 


10T=0 
DO 53 I=1,NRU 
1F (FL (13.£0.0.0) GOTO 53 
10T=1 
53 CONTINUE 
ІҒ (10Т.Е0.1) 6010 54 
WRITE (6,117) 
117 FORMAT (1X, ‘ххх RED FORCE IS ELIMINATED. END OF BATTLE. ') 


GOTO 66 

54 10T=0 
00 55 I=K,L 
ІҒІҒІ (1) .Е0.0.0) 6070 55 
10T=1 


55 CONTINUE 
ІҒ(10Т.Е0.1) 6010 65 
WRITE (6,118) 

118 FORMAT (1X, ‘xxx BLUE FORCE 15 ELIMINATED. END OF BATTLE. °) 
GOTO 66 


p 





551 6000 WRITE (6,119) 


552 119 FORMAT (1X, "ммм DISTANCE BETWEEN FORCES IS TOO CLOSE. END OF BATTLE 
553 17) 

554 GOTO 66 

555 65 1С=1С+1 

556 GOTO 67 

557 66 STOP 

558 END 
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SUBROUTINE LOS (XA, TA, THACA, TMICA, SIZEA, XB, 18. TMHACB, TMHICB, SI ZEB, 


-LRTOB,LBTOR, VISFRA, VISFRB) 


COMMON /HILLS/ XC (100) , YC (100) , PEAK (100) , Sx (100), SY (100) , RHG (190) 
COMMON /HILLS/ SCALE (100) , THGRHO (100) , THOSCL (100) , BASE 

COMMON /HILLS/ NHILLS 

COMMON /COVER/ CXC (150) ,CYC (150) ,CPEAK (150) ,CPXX (150) ,CPYY (150) 
COMMON /COVER/ CPXY (150) ,NCVELS 

COMMON /COUNTR/KH, KHN, KV. KN, KGRS, KELL, KINT 

COMMON /GRIO/ LST (10,10) , NHL (10, 10) , LISTH (50) , KHREP (100) , KTREP 
COMMON /GRID/ LSTC (10,10) , NC (10,10) ,LISTC (40Q) , KCREP (190) 
OIMENSION I6X(100),16Y (100) , IEL (100) ,CS1 (100) ,CS2 (100) 

DATA NGRIO/10/,GSIZE/1000./ 


С SUBROUTINE TO COMPUTE FRACTION VISIBLE FOR OBSERVER TARGET PAIRS 


VISFRA=1. 

VISFRB«1. 

XBRs XB- XR 

YBA=YB-YA 

IF ((XBR.EQ0.0.).AND. (1BA.£0.0.)) RETURN 
IF (SIZER+TMICA.LE.O.) GO TO 510 

IF (SIZEB*TMICB.LE.0.) GO TO 510 

1F (TMICA.LT.0.) VISFRA=1.0+TMICA/SIZEA 
IF (THICB.LT.O.) VISFRB=1.0+TMICB/SIZEB 
ZA=sTMACA + TMICA + SIZEA 

2ВеТИЯСВ + TMICB + SIZEB 

KTREPSKTREP*1 

2ВЯ-28-28 

XBASG=XBAxXBA 

YBASQ=TBA=TBA 

XYBA=XBAxTBA 

TWOXBAS2. xXBA 

TWOYBAs2. «YBA 


C COMPUTE 6810 SQUARES CROSSEO BY A TO B LINE 


95 

100 
110 
120 
125 
130 


140 


150 


NGRSO=0 

ІР (ХВЯ) 110, 95,100 
ХВАз0. 1 

156Х4-1 
XINC=GSIZE/XBA 

GO 19 120 

156X=1 
XINC=-GSIZE/XBA 

ІР (ТВА) 140, 125,130 
YBRsD.1 

1567+-1 
VINC=GSIZE/YBA 

G0 TO 150 

150141 
YINC=-GSIZE/YBA 

IXe te] F 1X (XB/GSIZE) 
1F (1X.GT.NGRIO) IX=NGRIO 


ol 





100 


IY*1*IFIX(YB/GSIZE) 
IF (IT. GT. NGRID). IY«sNGRID 
XNEXTzGSIZEx (FLORT (1X) *0. 5x CIS6X-1.)) 
YNEXTzGSIZEx (FLORT (1Y) *0. 5x CISGY-1.)) 
XSTEP= (XB-XNEXT) /XBA 
YSTEP= (YB- YNEXT) /YBA 
160 NGRSQ=NGRSQ+1 
IGX (NGRSQ) = IX 
ІСТ (NGRSQ) =1Y 
1F ((XSTEP.GT.1.) . AND. (YTSTEP.GT.1.)) GO TO 200 
IF (XSTEP-YSTEP) 170,180,190 
170 IXsIX %156Х 
XSTEP2XSTEP*XINC 
GO TO 180 
180 IXsIX*ISGX 
XSTEP=XSTEP*+XINC 
190 171 7+ 1567 
YSTEP=YSTEP+YINC 
GO TO 160 
200 KGRS=KGRS+NGRSQ 
C GRID SQUARE LIST NON COMPLETE IN IGX, IGY WITH NGRSQ ENTRIES 
С 
C NOW FIND WHICH COVER ELLIPSES TOUCH THE A TO B LINE, 
C CHECK ELEVATIONS AT S1 AND Se FOR EACH SUCH ELLIPSE 
МЕ: 520 
CHTMAX=0. 
IFINCVELS.EQ.0) GOTO 270 
DO 260 K=1,NGRSQ 
IXsIGX tK) 
1 та бтк) 
N=NC (IX, IT) 
IF(N.E0.0) 060 Т0 250 
LSsLSTC (IX, IY) 
LEND=LS+N- 1 
DO 250 L=LS,LEND 
KELL=KELL*+1 
ICsLISTC (L) 
IF (KCREP (IC) .EQ.KTREP) GO TO 250 
KCREP (IC) sKTREP 
RXsXR-CXC (1C) 
RY=YA-CYC (IC) 
PXX=CPXX (IC) 
PYY«CPYY (IC) 
PXY=CPXY (IC) 
RRsPXXxXBRSQO*PYT«YBRSQ*PXTXXTBR 
BB=PXXxTWOXBAxRX*PY YxTWOYBAXRY +PXY™ (RXXYBAtAY=XBR) 
CC»PXXxRXxRX*PYYx«RYxRY*PXYxRXxRY-1.0 
ARG=BB»*BB-4. OXAAXCC 
IF (RGRG.LE.0O.) GO TO 250 
$Q=SQRT (ARG) 
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101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
11Ч 
115 
116 
117 
118 
119 
120 


- 121. 


122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
14Ч 
145 
186 
147 
148 
149 
150 


C NOW 


210 


220 


230 


240 


250 
260 
C ALL 


C NOW 
270 


51е- (88:50) / (2. Ожйй) 

522 (S0-BB) / (2.О»йй) 

15 (51.0Е.1.) 00 10 250 

IF (S2.LE.0.) GO TO 250 
IF(SI.LE.0.) GO TO 510 

IF (S2.GE.1.) GO TO 510 

CHECK LOS RT S1 RND 52 
KINT=KINT+1 

CPKsCPERK (1C) 

XS2X8*S2xXBA 

YSsTYAtS2xYBA 

CRLL ELEV (XS, YS, HTS) 
HTSeHTS*CPK 

Z25228*52«2BR 

IF (LATOB.EQ.0) GO TO 210 
CALL KOVER (ZA, TMACB, SIZEB, 2B, $2,HTS, ZS, VISFRB) 
IF (VISFRB.LE.0.) GO TO 510 
IF (LBTGA.EQ.0) GO TO 220 
Se1.-S2 

CALL KOVER (2B,TMACA,SIZEA,ZA,S,HTS,ZS, VISFRA) 
IF (VISFRA.LE.0.) GO T8 510 
XS«XR*S1xXBR 

Ү5а(й»51хҮВВ 

CALL ELEV (XS, YS,HTS) 
HTS=HTS+CPK 

29228*S1«2BR 

ТЕ цатав. Е0.0) GO То 230 
CALL KOVER (ZA, TMACB, SIZEB, 28.51,HTS, 25, VISFRB) 
IF (VISFRB.LE.0.) GO TO 510 
IF (LBTOR.EQ.0) GO TO 240 
ба1.0-51 

CALL KOVER (2B,TMACA,SIZEA,ZA.S,HTS,ZS, VISFRA) 
IF (VISFRA.LE.O.) GO TO 510 
NELS=NELS?1 

IEL (NELS) «1C 

CS1 (NELS) =51 

CS2 (NELS) a52 

IF (CPK.GT.CHTMAX) CHTMAX=CPK 
CONTINUE 

CONTINUE 

ELLIPSES CHECKED 


START ON THE HILLS 

DO 600 Ks1,NGRSQ 

1Х«16Х (К) 

| Та 1 6Т (К) 

IF (NHL (IX, 17) .EQ.0) GO TA 600 
LSsLST (1X, IY) 
LENDSLS*NHL (IX, IY) -1 

09 500 L=LS,LEND 


65 





151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
139 
200 


IsLISTH (L) 
IF (KHREP (1) , EQ. KTREP) GÓ TO 500 
KHREP (1) «KTREP 


C PROCESSING FOR HILL I STARTS HERE 


Кнакн 1 


C COMPUTE W «ТОР OF THIS HILL ALONG O-T LINE 


280 


390 


СХаХВЯ/5Х (1) 

CY#YBA/SY (7) 

DXs (XR- XC (12) /SX (I1) 

DYs (YR-YC (1)) /SY (1) 

FQ=TWOSCL (1) * COXXDX*CYXxDY*RHO (1) * CXxQY *CY«DX) ) 
GO=SCALE (I) x (CCXx«CX*CYXxCY* TNHORHO (1) «CX«CY) 

IF (GQ .£Q.0.) GO TO 500 

N=-FQ/ (2.60) 

ТЕ (ABS (WN) .GT.5.) GO TO 500 

FSQ=FQxFQ 

EQ=SCALE (1) м (ОХм0х + ОТжОТ + ТИОВНО (1) xDXxDY) 
POWER=EQ-FSQ/ (4. GQ) 

IF (POWER .LT. -3.) GO 190 500 
HHN=PEAK (1) «EXP (PONER) 

KHW=KHW? | 

IF (HHN.LE.BASE) GO 10 500 

ZW2ZAtWNZBA 

IF((N.LT.0.).0R. (W.GT.1.)) GO TO 300 

IF (HHW .GE. ZW) GO TS 510 

CvHTH=0. 

IF (N9ELS.EQ.0) GO TO 300 

DO 280 M=1, NELS 

IF ((CS1 (M) .GE.N) .OR. (CS2 (M) .LE.N)) GO TO 280 
ICsIEL (M) 

IF (CVHTH. LT. CPERK (100) CVHTNSCPERK (1C) 
CONTINUE 

IF (CHHW*CVHTW) .GE.ZW) GO TO 510 

IF (HHN*CHTMAX.LT.AMINI (ZA-SIZEA, ZB-SIZEB)) GO TO S00 


C IF NE GET TO HERE THEN NEEO TO FINO LOWEST SIGHT LINE OVER HILL 
С WEWTON ITERATION A TO B GIVING VISFRB 


330 


IF 1.АТӨВ.Е0.0) GO TO 400 

Ку-ку? 1 Í 

Уан 

HHV=HHN 

NCT=0 

FY=FQxY 

ТНОбУ=2. мбОму 
ЕСМух?7А + ННУж (ТНОСУжу + РУ-1.) 

КМ=КМ+ 1 

FACTOR= (THOGVATHOGVU +2. (GO+THOGVxFQ)+FSOQ) 

DF CNVzHHVx«VxFACTOR 

IF (ABS (DFCNV) .LT.1.E-10) GO TO 350 

Y2V-FCNV/DFCNY 

IF (ABS (V) .GT.5.)GG0 TO 400 





201 
ede 
203 
208 
205 
206 
207 
208 
209 
210 
211 
212 
213 
elu 
215 
216 
217 
218 
213 
220 
221 
егг 
223 
224 
225 
226 
гг? 
228 
гез 
230 
231 
232 
233 
234 
235 
236 
237 
238 
239 
280 
241 
eue 
283 
гчч 
245 
246 
247 
288 
289 
250 


350 


380 
390 


FV=FOxV 

TWOGV22. xGQxV 

POWER » EQ*FV*GüxVxV 

IF (POWER .LT. -3.) GO TO 400 
ннугРЕЙВК (1) хЕХР (РОМЕН) 
Оннуаннух (Ей + ТНОСУ) 

Е | уа? А + ОННУжу 

IF (ABS (HHV-ELV) .LT.1.) GO TO 350 
NCT=NCT+1 

IF (NCT.LT.10) GO TO 330 

IF C(V.LT.0.).0R. (V.GT.1.)) GO TA YOO 
СУНТУ=0. 

IF (NELS.EQ.0) GO TO 390 

DO 380 M=1,NELS 

IF ((CS1 (M) .GE.V) .OR. (CS2 (mM) .LE.VIIGO TO 380 
ICsIEL (M) 

IF (CVHTV.LT. CPERK CIC) ) CVHTVsCPERK (1C) 
CONTINUE 

HTVsHHV*CVHTV 

ZV=ZA+V"ZBA 

CALL KOVER (ZA, TMACB, SIZEB,ZB,V,HTV, ZV, VISFRB) 
IF (VISFRB.LE.0.) GO TO 510 


C NEWTON ITERATION B TO A GIVING VISFRA 


400 


430 


450 


IF(LBTOR.EQ.O) GO Т0 500 

KVsKV*1 

Уан 

ҮМізу-1. 

ННУзННИ 

NCT=0 

Руа хү 

TWOGV=e.xGQxV 
FCNV=ZBtHHV™ ( (FO+TWOGV) xVM1-1.) 
KN=KN+ 1 

FACTOR= (THOGV=TWOGVt2. x (GO+ TNOGVxFQ) *FSQ) 
DFCNV2HHVxVM1 «FACTOR 

IF (ABS (DFCNV) .LT.1.E-10) GO TO 450 
VaV-FCNV/DF CNV 

IF (ABS (Y) .6T.5.)G60 TO 500 

ум! зу-1. 

FYsFQxV 

TNOGVa2.xGQxV 

POWER = EQtFVtGQxVxV 

IF (POWER .LT. -3.) GO TQ 500 
HHV=PEAK (1) xEXP (POWER) 
Оннуаннух (Е0+ ТНОСУ) 

ELVsZB*DHHVXxVMI 

IF (ABS (HHV-ELV) .LT.1.) GO TO 450 
NCT=NCT+1 

IF (NCT.LT.10) GO TO 430 

IF ((V.LT.0.).0R, IV.GT.1.)) GO TO 500 


02 





251 
252 
253 
25Ч 
255 
256 
257 
258 
253 
260 
261 
262 
263 
264 
265 
266 
267 
268 
269 


480 
430 


500 
600 


510 


CVHTV=0. 

IF (NELS.EQ.0) GO TO 490 

DO 480 M=1,NELS 

ІЕ ((С91 (М). СЕ. У) . ОА. (С52 (М) .LE.V))GO TO 480 
ICsIEL (M) 

IF (CVHTV. LT. CPERK (1C) ). CVHTVsCPERK (1C) 

CONT INUE 

НТУ< ННУ + СУНТУ 

ZV=ZAtVxZBAR 

S2-VM1 

CALL KOVER (ZB, TMACA,SIZER,ZA,S,HTV,ZV, VISFRA) 
IF (VISFRA.LE.O.) GO TO 510 

CONTINUE 

CONTINUE 

RETURN 

VISFRA=0. 

VISFRB=0. 

RETURN 

END 
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о юг ». 


uu 
45 
46 
47 
ча 
48 


u? 


17 


50 


37 


27 


60 


65 


SUBROUTINE SETUP 


COMMON /HILLS/ XC (100), YC (100) , PEAK (100) ,SX (100) ,S¥ (100) , RHO (100) 
COMMON /HILLS/ SCALE (100) , THORHG (100) , TROSCL (100) , BASE 
COMMON /HILLS/ NHILLS 

COMMON /COVER/ CXC (150) ,CTC (150) ,CPEAK (150) ,CPXX (150) ,CPYY (150) 
COMMON /COVER/ CPXY (150) , NCVELS 

COMMON /COUNTR/KH,KHH,KV,KN, KGRS, KELL, KINT 

COMMON /GRIO/ LST (10,10), NHL (10,10) ,LISTH (450) , KHREP (100) , KTREP 
COMMON /GRIO/ LSTC(10,10),NC (10,10) , LISTC (400) , KCREP (150) 
RERD (8,7) NHILLS 

READ (8,47) BASE 

FORMAT (F 10.4) 

FORMAT (16) 

FORMAT (SF8.3,F6.4) 

08 SO Is1,NHILLS 

READ 18,17) ХС {1) ‚ТС (1) PEAK (1) ,SX (1) ,SY (1) ,RHO 1]) 
CONTINUE 

READ (8, 37) LST 

READ (8, 37) NHL 

READ (8, 7) NHTOT 

READ (8,37) (LISTH (1) ,1=1,NHTOT) 

FORMAT (1615) 

READ (8, 7) NCVELS 

IF (NCVELS.EQ.O) GO TO 65 

00 60 J=1,NCVELS 

READ 18, 27) CXC (1) ,CYC 11) ,CPEAK (1) ,CPXX (1) ,CPYY (1) ,CPXY (1) 
FORMAT (3F10.4,3E 13.77) 

KCREP (1) =-214 7483600 

CONTINUE 

RERD (8, 37) L3TC 

RERO (8, 37) NC 

RERD (8, 7) NCTOT 

RERO (8, 37) tLISTC(1), 131, NCTOT) 

DG 100 I*1,NHILLS 

SX (1) =SX (1) «1.625 

SY (1) =SY (1) “1.625 

ХС (1) = (ХС (1) -500.) „100. 

ҮС (12 а (ҮС (1) -930.) «100. 

TWORHO (7) =2. xRHG (1) 

SCALE (1) =-1./ (2. (1. -RHO (1) x2) ) 

TWOSCL (1) =2. SCALE (1) 

KHREP (1) =-214 7483600 


C ALL VALUES NOW IN METERS ON O -- 10,000 GRID 


100 


CONTINUE 
KTREP=-2147483600 
КНа0 

KHW=0 

КҮз0 

KN=0 
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51 
52 
53 
5Ч 
55 


KGRS=0 
КЕ! з0 
КТ МТ 0 
RETURN 
END 
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SUBROUTINE ROUTE 


COMMON /GRP3/ NBU,NRU,FL (6) , FO (6) , NOI (32 ,X1C (3,200) ,YIC (3,200), 
1101R (3,200) ,AVSP, 15Р0 
1, IUSTRT (69) , I1 (B) , LOST (6, 8) , VISFRR, VISFRB, SIZETK, 
1S12€ TH, NT (6) , NF (8) , SRF, DISMAX, 
INLOSC (6,6) , VISFR (8,8) , RHINTK,RMXTK, RMINTH, RAXTH, OP, TONFR, TNKFR, 
IPTT (3, 3) , RF, PR (6,8) , RPOR (8,8) , LOR (8, 8) , NN (B) , OFL (8) , POL (6) 
DIMENSION XLOC (3,20), YLOC (3,20) ,N (3) 
ТЕ (15Р0.Е0.Ч) DST=80.463 
IF (1SPD.EO.3) DST=67.053 
ТЕ (15Р0.Е0.2) OST=53.643 
IF (15Р0.Е0.1) 057+40.232 
DO 300 1=1,NRU 
READ (5,15) N(1) 
15 FORMAT (12) 
NL=N (1) +1 
DO 200 IN=2,NL 
READ (5,201) XLOCS, YLOCS 
201 FORMAT (F6.1,1X,F6.1) 
XLOC (1, 1N) =XLOCS 
YL£C (1, 1N) =YLOCS 
200 CONTINUE 
XL OC (1,1) =XIC (1, 1) 
YLOC(1,1)=Y1C (1, 1) 
1DIR (1,1) =0 
NL=N (]) 
NUM=2 
06 305 J=1,NL 
XLaXLOC (1, J*1) -XLOC (1, J) 
YL#YLOC (1, J* 1) - TLOC (1, J) 
DISTsSQRT (XL«x2* YL xx2) 
Y=ABS (YL) 
Z=Y/XL 
ANGL=ATAN (2) 
ОЕбаАамбі «57. 2958 
IF (J.EQ.ı) GO TG 320 
XLN= (DST-EXTRA) «COS (ANGL) 
DISTs (DIST*EXTRR) -DST 
YLNs (DST-EXTRR) «SIN (ANGL) 
X1C (I, NUM) =X1C (1, NUM-1) *XLN* XLE 
IF (YL.GT.0.) GO TO 325 
YLNz-YLN 
325 YIC (1,NUM) =Y1C (1, NUM-1) * YLN* YLE 
IF (YL.GT.0.) 60 TO 340 
1018 (1, NUM) =-1F 1X (DEG) 
00 ТӨ 341 
340 1018 (1, МУМ) =1Е1Х (0Еб) 
3401 NUMzNUM*1 


a? 





51 
52 
53 
SY 
55 


57 
58 
59 
60 
61 
62 
63 
64 
65 


67 
68 
69 
70 
71 
72 
73 


320 


310 


342 
343 


315 


305 
300 


XLN=DSTx*COS (ANGL) 
YLN=OSTSIN (ANGL) 

IF (YL.GT.0.) GO TO 310 
YLN2- YLN 

IF (DIST.LT.DST) GO TO 315 
X1C CI, NUM) »XIC (I, NUM- 1) *XLN 
YIC (I, NUM) 3YIC (1, NUM- 1) *YLN 
IF (YL.GT.0O.) 008 TO 342 
IDIR (1, NUM) »- IFIX (DEG) 

GO TO 343 

IDIR (I, NUM) 31F1X (DEG) 
DISTs01IST-DST 

NUM=NUM+ 1 

GO TO 310 

EXTRA=DIST 

XLESEXTRAxCOS (ANGL) 
YLE=EXTRASIN (ANGL) 
IF(YL.CT.0.) 60 TO 305 
YLE=-TLE 

CONTINUE 

CONTINUE 

RETURN 

ENO 
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SUBROUTINE LRMOR (1, J, PCTVIS, OETRRT,P 


K 
SUBROUTINE TG COMPUTE DETECTION 0008 OF TARGET J 
BY OBSERVER 1 GIVEN THE VISIBLE FRACTION (PCTVIS). 


COMMON /GRP1/ IPROIA (6) , ISECHD (6) ,MVTOIR (6) , X (6) , Y (6) , SPD (6) 
TCFACT=1.0 

ZEROL=0.00001 

PAIs22.0/7.0 

02 (15ЕСНО (1) »РАТ / 180. 0) /2.0 

BBB= (1.0/7 (2.0% (SIN (0) -DxCOS (D)))) 

IF (ABS (BBB) .LT.ZERÓL) BBB=0.0 

RRR» (- BBB) «COS (0) 

IF (ABS (ARA) .LT.ZEROL) ARA=0.0 
OTANG=ATANE ( (Y (3) -1 (19), (X (3) -X (133) 
POsIPROIR (1) <PA1/180.0 

IF ((PDxOTANG) .GE.0.0) GOTO 1 

IF (PD.LT.0.0) GOTO 9 
ANGLE=2xPA1+0TANG-PD 

GOTO 10 

RNGLEs2xPRI*PO-OTRNG 

IF (ANGLE. GT. PRI) ANGLE=2xPA1-ANGLE 
GOTO 2 

ANGLE=ABS (PO-OTANG) 


2 JF (RNGLE.GT.D) GOTO 3 


6 


DUP=PD+D 

DLOWN=PD-D 

RNGLF TSOTRNG* (15. Ox«PRI/180.) 

IF (ANGLFT. GT. OUP) ANGLFT=0UP 
RNGLRTsOTRNG- (15. xPRI/180.) 

IF (ANGLRT.LT. OLON) ANGLAT=OLOW 

PK2BBBx (SIN (ANGLF T) - SIN (ANGLRT) ) *RRR* (ANGLF T-RNGLRT) 
IF (PK.LT.0.0) GOTO 3 

IF (PK.GT.1.0) GOTO 5 

6070 8 

РКа0.0 

OETRAT=0.0 

0010 6 

PK=1.0 

RANGE=SQRT ((X (9) -X (1) ) жже (7 (.)) -У (1)) мма) 
АА =0. 001 «ВАМСЕ/РСТУ15 

TORNG=ATAN2 ( (Y (1) -1 (33), (X (1) -X ())) 
ЯО=МУТОТЯ (.}) »РАТ/ 180. 0 
HORVEL=ABS (SPD (J) *SIN (TÓRNG-RO) ) 
DENOM=1.453+TCFACTA (0.5978+2. 188 (RAxx2) -0. 5038xHORVEL) 
IF (DENOM.LE.ZEROL) OENOM=ZEROL 
DETRAT=0.003+1.088/0ENOM 
DETRAT=OETRAT=PK 

RETURN 
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C FUNCTION TO COMPUTE TERRAIN ELEVATION FOR GIVEN X, 


100 
1S0 


SUBROUTINE ELEV (X,Y, TMAC) 


COMMON /HILLS/ XC (100) , YC (100) , PERK (100) , SX (100) , ST (100) , RHQ (100) 
COMMON /HILLS/ SCALE (100) , THORHO (100) , TWOSCL (100) , BASE 


COMMON /HILLS/ NHILLS 


COMMON /GRIO/ LST (10,10) , NHL (10, 10) , LISTH (4S0) , KHREP (100) , KTREP 
COMMON /GRID/ L3TC (10,10) , NC (10, 10) ,LISTC (400) ,KCREP (150) 


DATA NGRID/10/,GS1ZE/1000. / 


ZMAX=BASE 
IX21*1F IX (X/GSIZE) 

IF (IX. GT. NGRID) IXsNGRID 
IY=1+1F1X (1/GSIZE) 

IF (IY.GT.NGRID) IY=NGRID 
IF (NHL (1X,1Y).£Q.0) GQ TO 150 
LS=LST (1X,1Y) 
LENDsLS*NHL (IX, 11) -1 

DO 100 L=LS,LEND 

T=LISTH (L) 

QX» (Х- ХС (1) ) /SX (1) 
0х59=@хХмах 

IF (QXSQ .GE.9.) GO TO 100 
QY» (Y- YC(12) /SY t1) 
QYSQzQY«QY 

IF (QYSQ .GE. 9.) GO TO 100 
QXYs TNORHÓ (1) xQXxQTYT 
FACTOR=SCALE (1) x (QXSQ*QYSQ-*QXT) 
IF (FACTOR.LT.-3.) GO TQ 100 
HT=PEAK (1) «EXP (FACTOR) 

IF (HT.LE.ZMAX) GO TO 100 
ZMAX = HT 

CONTINUE 

TMAC=ZMAX 

RETURN 

END 
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Y COORDINATES, 





1 SUBROUTINE STGCH (I, RANGE,A) 

2 E 

3 С SUBROUTINE TO COMPUTE STOCHASTIC ATTRITION COEFFICIENT 

4 C 

5 COMMON /GRP6/ ALPHA (6) 

6 COMMON /GRP3/ NBU, NRU,FL (6) , FO (6) , NO1 (31, XIC (3,200), TIC (3,200), 
? 1IDIR (3, 200) , RVSP, ISPD 

8 1. JUSTAT (6) ,11 (6) ,LOST (6,6) ,VISFRA,VISFAB,SIZETK, 

9 1SIZETH,NT (8) ,NF (6) ,SRF,DISMAX, 

10 INLOSC (6,61 , VISFR (6,6) ,RMINTK, RMXTK, RMINTN, RMXTN, OP, TONFR, TNKFR, 
11 ІРТТ(3,3),АҒ,РОЯ (6,6), АРОЯ (6,6),108 (6,6! ,М8Я (6) ,ӨҒІ. (6) , POL (6) 

12 ІЕ (1.Е0.22 60 TO 5003 

13 As ALPHA (1) « ((1.0-RANGE/RMXTH) x2) 

14 GO Ta 5004 

15 S003 A=ALPHA (1) = ((1.0-RANGE/RMXTKI 2) 

16 5004 RETURN 

17 ENO 


e 





Фо O <J O) U1 YN ~ 


SUBROUTINE ETK(I,RRNGE, T) 
SUBROUTINE TO COMPUTE EXPECTED TIME TO KILL A TARGET. 


COMMON /GRP2/ TR(2),T1 (2), TH (2), TM (2), TF1 (2) , TF2 (2), TF3(21, 


IP (2,60) , PHH (2,60 , PHA (2,68) , PKH (2,6) , TF (2) 


IF (I.EQ.2) GOTO S 

TF (1) #TF1 (1) 

GOTO 6 

IF (RANGE.GT.1000.03 6010 7 

TF (I) *TF 1 (1) - (TF 1 (12 « (2000. 0- RANGE) /1000.0) 

GOTO 6 

IF (RANGE.GT.2000.0) GOTO 8 

TF (1) =TF2 (1) - ((TF2 (1) -TF1 (13) « (2000. 0-RANGE) 1000.0) 
сото 6 

TF (1) 21F3 (1) - ((TF3 (1) -TF2 (1) ) x (3000.0-RANGE) /1000. 0) 
Js (RANGE +250.0) /500.0 

IF (J. GT.6) J56 

TsTR (CI) *T1 (1) - THCI) * CCEH COH) * TF (02) /PKH C, 02) * CCEM CD * TF (2) 7 


IPHM C1, J2 ) x CCI.O-PHH C1, ? /PKH C1, 3) +РНН (1,4) -P (1, J3) 


RETURN 
END 
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won ook wh = 


21 


19 


SUBROUTINE SORT (1,m) 

COMMON /GRPS/ LOT (6,6) ¿ROT (6, 6) 
DO 19 J=1,M 

IF (ROT (1, 0 .GE.ROT (1,J)) GOTO 19 
R=ROT (1, J) 

NN=LOT (I, J} 

ВОТ (1,.) «В0Т(1,М) 
LOT(1,J)=LOT(I,M) 

ROT(1,M)=R 

LOT (1,M) =NN 

CONTINUE 

RETURN 

END 
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wan num E Wf = 


2000 
2050 


SUBROUTINE KOVER (20, THACT, SIZET, 21,5, HTS,25, VISFRT) 


ЇГ (5.Е4.0.) 68 10 2000 

IF (HTS.GE.ZS) GO TO 2050 
НЕХТг20% (НТ5-10) /5 
EVIST#AMAX1 (HEXT, TMACT) 

IF (EVIST.GE.ZT) GO TO 2050 
IF (EVIST. LE. ZT-S5IZED RETURN 
VISs (ZT-EVIST) /SIZET 

IF (VIS.LT. VISFRT) VISFRTsVIS 
RETURN 

IF (HTS.LT. Z0) RETURN 
VISFRT=0.0 

RETURN 

END 
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APPENDIX С 


Definition of Variables in Computer Program 


ALPHA(I) 


APOA(I,J) 


AVSP 
BREAK 


DISMAX 


DIST 


DST 


22001) 
FO(T) 
LALT 


ІС 


IDIR(I,J) 
ШШ (Т) 
ІМОУЕ 


IPRDIR(I) 
TRTE 


IS 
ISECWD(T) 


Шог әрігісісп-гасе coefficient for 
Seoehastıc atverıitiıon module, 


Ihe average proportion of the 1 n attacker 
eunt a al locat nd 10 ТЇїгсө 0п unit i. 


Average speed of moving attacking units. 


Breakpoint distance between attackers and 
defenders. 


faximum distance allowed between attacking 
units before the leading unit is delayed. 


The straight-line distance between two 
movement nodes inputed by the user. 


The distance to be moved each time step 
by attacking units. . .— 


Force level of unit i. 
Виета тосе ауе of unit le 


Denotes whether the user desires alternate 
defensive vositions or not, 


ОПП с number Of time units a defender 
has been moving. 


th routes 


Direction of j"? interval in i 
limecnval si ndex for unit i; 


Number of time units a defender is allowed 
for moving to an alternate position, 


Primary direction of movement for unit i. 


Denotes whether user wants to input routes 
От по» 


Киоссе сеа пед Гог stochastic attrition. 


ГІНЕН Search Sector for mit il. 


20 


Ни 








ESPD 


ІТЕМ 


[TIME 
РТТ 


Imam (I) 
O T.J) 
гозер 


I (T.J) 
MVTDIR(I) 
N(I) 


NA(I) 
NBU 
NF(I) 


ОБС (Т.Ј) 
HOI(I) 
NRU 


ШЕГІ) 
OF L(I) 


RI. LJ) 


i T.J) 


Input variable to denote user's desired 
sveed for attackers movement. 


Input variable denoting number of time 
steps allowed for defender's move. 


Ситлештттъй те, Seconds, of battle, 


Input variable denoting whether attrition 
ШЕ: be SvOChastic or deterministic. 


Current Statúüs of unit 1% 


h 


The number of the jt attacker or unit i.e 


Denotes whether line-of-sight exists 
between unit i and j or not. 


th 


The number of the j ее. 


Movement direction of unit i. 


Number of nodes inputed by user for 
011621. 


oDer eo attackers. of- unit i- 
Number of blue units. 


Number o time units unit 1 is allowed to 
fire sat the same locations. 


Number of continous time steps that 
line-of-signt does not exist between unit 
Far Unie je 


de 


Number of intervals in the ЕНЕ route, 


Number of Red Units. 
Number ot targets of unit i. 


Force level of unit i during previous 
time step. 


Probability of 19" 


in range band j. 


ӘЛПЕТ by unit- i 


Probability of a hit following a hit 
by unit i in range band j. 
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PHL(I,J) 


ESSCT,J) 


PH 


POA(I,J) 


POL(I) 


mit, J) 


RANGis 


RMINTK 


ETA 


КМХТК 


RHXTW 


ESTCI,J) 
ВИК 
SIZET\I 
TA(K) 


ee CK) 


BECK) 


TES(K) 


Probability of a hit followins a miss 
œ unit 1 іп range band 1. 


вора ел о а kill given a hit by 
Una range band jJ. 


ШИЕ 220101 027 time a moving unit 
is searching for targets. 


The proportion of the j*" attacker of 
ит а осасеа to fire on unit i. 


Percento unie i lost during the 
current time sten. 


ео о A ving firepower 
allocated to the it. target if there 
are j targets available. 


Current minimum distance between attackers 
and defenderse 


Kinimum effective ranse for attacking 
weapon system. 


Minimum effective range for defending 
weapon system. 


Maximum effective range for attacking 
weanon system. 


Maximum effective range for defending 
weapon system. 


The range of the n 


Coe COl unit. Ls 
size of attacking vehicle, 
Size of defending vehicle. 


tame to acguire a target for us weapon 
System type (k = 1,2). 


о ТЕКЕ 50 -1000т tor кей weapon 
system type (k = 1,2). 


шил СН ІСЕР Lo 2000m Tor js weapon 
system type (k - 1,2). 


Mame sot СЕБЕ ko »000m for „en weapon 
Биеп буре (К = 1,2). 


Да 





ТН(К) 


ВОСК) 


ТМ(К) 


TNKFR 
ТОШЕВ 
COLI) 


VISFR(I,J) 
VISFRA 
VISFRB 
хаттаар 
ГЕЛ). АСТ) 
Pe T.J) 
СО, 
PSY T, 


ВОС (т, 
ВОС (Т, Ј) 


Time to fire a round following a hit 
for weapon system type k (k = 1,2). 


ШЕСЕ со.  г  О0па а1Ссеёег target 
has been acquired for weapon system 
001027. 


meto tire a round following a miss 
for weapon system type k (k = 1,2). 


Firing rate for attacking weapon system. 
Firing rate for defending weapon system. 


Total vercentage of lost since battle 
began for unit i. 


PoS pact Ontor mnit i seen by unit- 1. 
Fraction of unit A as seen by unit &. 

Шоо ст unit Bas seen by unit A. 

Coordinates of unit i. 


Coordinates of alternate position for 
defender i. 


Coordinates of the jun interval endpoint 
оле route for unit il; 


Distance added to previous interval 
endpoint for vehicle to move DST during 
a time stepe 


Coordinates of the ule node inputed by 
Pre оо tor che route of unit i. 
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